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ABSTRACT 

A compact heat exchanger testing facility was cbnstructed at the 
United States Naval Postgraduate School. The testing facility is used 
to collect heat transfer and friction factor data by the steady state 
steam to air technique. Heat exchangers with frontal dimensions of up 
to 12 inches by 12 inches can be tested. The maximum pressure drop which 
can be produced across the air side of a heat exchanger is 2.5 psi. The 
testing facility was designed to obtain data for exchanger air side 
Reynolds numbers from 500 to 10,000. A description of the apparatus, the 


experimental method and equations and an evaluation are given. 
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NOMENCLATURE 


Letter Symbols 


Air side total heat transfer area, £t- 


Air side free flow area, ft? 


Air side fin area, ft 
Air side frontal area, ft? 

Steam side fin area, ft2 

Steam side total heat transfer area, ft2 
Air side wall area, ft2 

Steam side wall area, ft 
Coefficient of discharge for a metering orifice 


Jet contraction ratio 


Coefficient of discharge of a fluid orifice for a given value 
of 8B 7 


Specific heat of air at constant pressure, BTIU/(1b, °F) 
Diameter of fluid metering orifice, in. 
Air orifice duct diameter, in. 


Hydraulic diameter of duct immediately downstream of test 
core, in. 


Error in heat balance between steam side and air side, percent 
Fanning friction factor in test core 

Friction factor in duct immediately downstream of test core 
Friction factor of a passage in test core 

Velocity of approach factor for the fluid metering orifice 
Thermal expansion factor for the fluid metering orifice 


Proportionality factor in Newton's Second Law, acm 32.174 
lb, ft/(1b¢ sec¢) 


Test core air mass velocity, 1lbs,,/ (hr ft2) 


VTi 





Unit conductance for thermal convection heat transfer 
(film coefficient), BTU/(hr ft2 °F) 


Enthalpy of the condensate leaving the test core, BTU/1b,, 
Enthalpy of steam, BTU/1b, 

Humidity ratio for air, 1b, water/lb, dry air 

Unit thermal conductivity, BTU/ (hr ft2 OF/ft) 


Contraction coefficient 


Velocity distribution coefficient 


Expansion coefficient 
Fin length equal to one half: of total fin height 


Length of duct from downstream face of test core to downstream 
pressure tap, in. 


Wall thickness between steam and air side of test core, in. 
Test core flow length, ft. 

V2 h/ké 

Mass rate of air flow, 1b,/hr. 

Pressure, in. Hj0, in. Hg, lb¢/in2, lb, /ft? 

Barometric pressure, in. Hg 

Heat transfer rate, BTU/hr 

Gas constant for air, 53.3 ft lb¢/(1b, °R) 

Hydraulic radius, (A.L/A,) , ft, (4 rp-hydraulic diameter) 
Mass rate of steam flow, 1b ,/hr 

Temperature, OF 

Absolute temperature, 9R, OK 

Unit overall thermal conductance, BTU/(hr °F ft2) 
Velocity, ft/sec 

Specific volume, £t?/1b_, 


Mass rate of condensate from test core, 1b} ,/hr 


vil 





Ratio of pressure upstream of the fluid metering orifice to 
pressure differential across the orifice 


Humidity correction to the specific heat of air 
Humidity correction to the density of air 


Net expansion factor for a square-edged metering orifice 


Greek Letter Symbols 


B 


A 


Meangdensity or air im test core for hot core test, Ibm / ie 


Ratio of fluid metering orifice diameter to duct diameter 
Denotes difference 

Fin thickness, in. 

Fin temperature effectiveness 

Total surface temperature effectiveness for air side 
Total surface temperature effectiveness for steam side 
Ratio of free-flow to frontal area, A ./AgR 

Dynamic viscosity, lbm/(hr ft) 

Density, lb +)/ft? 


Average density of air in test core for isothermal test, 
ibe ee 
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Dimensionless Groupings 


NRe 
Not 


Npy 
j 


NTU 


Reynolds number, (4r,G/4) 
Stanton number, (h/Gcy) 
Prandtl number, (Arc, /k) 
Colburn-j (Nop Np, 2/3) 


Number of heat transfer units for an exchanger (A, U/¢ Cp) 
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Subscripts 
l = 


2 = 


so - 
sl - 


s2 - 


Air upstream of core 
Air downstream of core 
Air metering orifice 
Air side of test core 
Core 

Fin 

Air film in test core 
Steam side of test core 
Steam metering orifice 
Steam in top header 
Steam in bottom header 
Wall between air and steam sides of test core 
Square fin 

Triangular fin 


Circular tube 


Force and Mass Units 


Ib. a 


lb¢ - 


Denotes pounds mass in distinction to 


Denoting pounds force 


ix 





ie Introduction 

Steady state tests transferring heat from condensing steam to 
flowing air are used to determine experimentally the heat transfer 
characteristics for compact heat exchangers. The friction factor is 
determined from the air pressure drop across the heat exchanger. Steam 
and air are used as the working fluids in the tests because they are 
readily available and produce effects which are easily analyzed. The 
heat transfer characteristics are given by Colburn-j heat transfer 
modulus evaluated from the steady state steam to air test (referred to 
here-in as hot core test). The friction factor is evaluated from both 
Enemnot cOre test and the tsothermal test (referred to as cold core test). 
The Colburn-j heat transfer modulus and friction factor are normally 
presented graphically as functions of the heat exchanger Reynolds number. 
Development of these three parameters is given in the theory section. 

In addition to determining the Colburn-j and friction factor, a 
heat balance is performed on the exchanger to check the consistency of 
the experimental data. 

The testing facility is similar in design to one built at Stanford 
University in 1947. The Stanford facility was built under a United 
states Navy contract to evaluate finned surfaces for use in compact heat 
exchangers and is described by Kays in reference [3]. 

Compact heat exchangers which are suitable for steady state steam 


to air testing have plate-fin geometry as shown beloy. 








The steady state steam to air testing facility was designed and 
built to accommodate a heat exchanger with frontal dimensions of up to 
12 inches by 12 inches. The facility is capable of producing air side 
Reynolds numbers from 500 to 10,000. An air side pressure drop of 2.5 
psi can be obtained. 

The testing facility is described and the experimental method and 
equations for data reduction are presented. The data reduction equations 
are conveniently arranged for programing on a digital computer. A test 
of a heat exchanger of known characteristics is presented as an evalua- 
tion of the testing facility. The operating procedures and the computer 
program for data reduction with accompanying instructions are given in 


the appendix. 





2. Theory 

Heat Transfer Characteristics. The heat transfer characteristics 
of a compact heat exchanger are best represented by dimensionless group- 
ings. This generalizes the results so that they may be used for working 
fluids other than steam and air. The succeeding development is that 
commonly used in the field; see Kays and London [4]. The Colburn-j heat 
transfer modulus is selected for this purpose as it combines heat transfer 
characteristics and fluid properties. The Colburn-j modulus is defined as: 
/ 2/3 
je Nop Nee’ = coe 

Cp kK 
where the unit conductance for convective heat transfer, h, must be deter- 
mined experimentally. The use of a steam to air system is one of the 
most direct methods of experimentally determining h. 

The unit conductance for convective heat transfer, h, for the air 
side, air being the fluid of interest, is determined from energy balance 
considerations within the axchanger. Low pressure slightly superheated 
steam, about 6 psig with 5 degrees superheat, enters the top of the heat 
exchanger steam side. The steam reaches its saturation temperature in a 
very short distance after entering the exchanger. Constant temperature 
is maintained as the steam condenses for the remainder of travel. The 
air is heated continuously as it passes through the exchanger. The 
sketch below illustrates the temperature variations in the working fluids 


as a function of the heat transfer area. 
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An energy balance for the differential area dA yields the equation 


for the heat transfer rate as: 


dq = hc, dTair = UdA AT 
(heat transfer (heat transfer rate (heat transfer rate 
rate) to the amr) from steam to air) 


Rearrangement gives: 


= dTair = UMA. 
(Tsteam-Tair) fn =D 


For small over-all changes in fluid properties and a relatively uniform flow 
distribution in the heat exchanger, the over-all unit conductance, U, 
may be assumed to be constant. Integration over the total area then 
yields: 

= Chama) = pte 
The dimensionless grouping on the right side of the equation is defined 
as the Number of heat Transfer Units, NIU. Thus the equation may be 


written as: 


gn ( Tsteam = Tair out = — NTU 
Tsteame] fair in 


Or 


( Tsteam - Tair out ) = oa 


Tsteam - Tair in 
The over-all unit conductance based upon the total air side heat 
transfer area must include a factor to account for the temperature 
variation in the fin areas. The over-all temperature effectiveness, n, 


may be written as: ie 
) = a — (1 -n¢) 
A 





where A¢ is the fin area andN Ff is the fin temperature effectiveness. 
The effectiveness for a fin of constant conduction cross section as 


given by Kays and London [4] is: 


ed tanh m& 
if mL 
2h : : 
where nm = kg for thin sheet fins 
2 = one half fin height for wall to wall fins 


Summing the resistances to heat flow then yields: 


1 = Ag Ag Lo + 1 

U ny Ag hy A, Ky No bh 
(over-all (steam side (wall (air side 
resistance) resistance) resistance) resistance) 


The convective coefficient for the air side, h, may be solved for giving: 





A value of h, equal to 2000 BTU/ (hr ft? °F) is assumed. A fifty 
percent error in this assumption results in approximately two percent 
error in the value of h. Accordingly the heat transfer coefficient for 
the air side, h, may be evaluated by experimentally measuring the fluid 
temperatures and the air mass rate of flow. 

Friction Factor. The flow of a fluid through an exchanger results 
in a pressure differential caused by the flow contraction as the entrance, 


the fluid friction in the exchanger and the flow expansion at the exchanger 





exit. These pressure effects are graphically represented below. 
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To define a meaningful heat exchanger friction factor, the entrance 
and exit pressure effects must be elimated from the total pressure 
differential. The friction factor will then be based solely upon the 
mean wall shear. The problem of entrance and exit pressure effects has 
been investigated both analytically and experimentally by Kays.[2]. 

The effect of pressure loss in the downstream duct is included in 
the total pressure drop in order that point 2 may be located downstream 
sufficiently to assure complete pressure recovery. The total pressure 


differential from point 1 to point 2 is then: 
Protal mele +OAPL - APH + APGuct 


The expression for the pressure effect due to flow contraction at 


the entrance is given by Kays [2] is: 


App ve (i- o*) 4 kevs 


¢ 


Was Age. 2c 





and for the flow expansion at the exit is: 


Ap, Ve (i- 7) _ keV 
Pb = qe aoc 


with g- denoting the exchanger free flow to frontal area ratio. K, dnd 
Ke are contraction and expansion coefficients which are functions of 
exchanger geometry and Reynolds number. Noting that the velocity may be 
related to the mass velocity and that the density change from section one 
to a and b to two will be very small compared to the total density change, 
the above equations yield: 


wed 


nae 2£4c Py 


ae c-2 + Ke) 


and 


——. ) 
App=zog (!- 02 — Ke) 
9c 8 
The pressure drop in the exchanger may be found by writing Bernoulli's 
equation in the differential form to include the Fanning friction term. 


Using the mass velocity, G, the equation as given by Kays [3] is: 


2 2 
cit— ac a GS {dx 
- ee 2 © 2gc@° (Th } 


~  - 
Zs 
Gc OC] Tb Ce. 
By defining the mean specific volume in the exchanger as: 
lL 
vm = (t) = 41 | dx 


and integrating from pn. QO to Ps L, the exchanger pressure differential 





Or 


is given by: 


é ec, fe Cs 
AP. = nm 2S =f eee | a 
c ede we ae = Th e 





The mean specific volume is found from its definition by replacing it 


with its value as given in the equation of state for a perfect gas: 


oe? = ge 
Vv RT 
and expressing temperature and pressure as functions of x. As the 


pressure variation is small compared to the temperature variation, an 


arithmetic average value of P may be used: 


The expression for temperature as a function of x may be written from 


the temperature-NTU relation as: 


f= to (Ts - Th) e 
where a is the steam temperature. Substitution of the expressions for 
temperature and pressure into the equation for the mean specific volume 


ields: 
a IL NTU ,. 
| 


2R 7 
Ven = = L(P3+P,) Tg - (1;-Ts} @ L dx 
O 


Integration and using the equation of state for a perfect gas to simplify, 
and assuming no temperature change from one to a and b to two the follow- 


ing equation obtains: 


iz 2 hae (T. — 7) 
PC, T3@s3 t+ Tblb NTU 


An alternate expression may be written by substituting the NTU-tempera- 


ture relation into the above equation. This expression is: 


g 2k Te — 
+ = £3 (T; —LMTD) 





with LMITD being the Log Mean Temperature Difference. 
The pressure loss in the downstream duct may be evaluated as: 
APyuct = 4 ie Vg 
Rd Dg 2c 

where the density, pj, the friction PACtOM@me eeche Length of duct 
from section b to 2, &q4, the hydraulic diameter, Dy, and the velocity, 
Va, are all quantities associated with the downstream duct. The velo- 
city in the downstream duct may be related to quantities appearing in 
the previous pressure loss equations. If the downstream duct has the 


same dimensions as the exchanger and small density variations are 


neglected, the continuity equation gives 


va (ENaC) 
Vy = a 








or 


Thence 
AP fosy G2 g-2 
= 4/——|} ——— 
duct | Dd 2dc Pr 
Again assuming that the density changes between section 1 to a and 


b to 2 are negligible, the a pressure drops are summed yielding: 


AP otal =s—— (joe at 2G al {i = &2 
C 





Manipulation of the equation gives: 





KO. ye ee oe 
total ~ ag. | Th &m 


f 
. (p+ ket £2 (1 +a) 


€2 





and solving the equation for the friction factor we obtain: 


7, Th &m Protal Ke cael) 
- 6¢/2q¢ eC; 
sd 
_ tek r o2(i+4 eH woe 


Co 


Therefore the friction factor, f, for a heat exchanger of given geometry 
may be found by experimentally determining the mass rate of air flow, 
the entrance and exit air pressures and temperatures and the steam 
temperature. 

Reynolds Number. For uniformity of reporting test results, the 
Reynolds number is defined in the same manner as given by Kays and 
London.{4]. This definition embodies the hydraulic diameter as the 
characteristic length and the velocity based on the minimum free flow 


area. Thus: x 
Ne, = Onve = Arnm® ~ The 
J 4 Ae y 
To evaluate the foregoing characteristics, the test apparatus 
must be capable of accurately measuring the air state upstream and 


downstream of the exchanger, the steam state within the exchanger and 


the mass rate of air flow. 


3. Description of the Test Apparatus 

General. Air is induced through one side of the compact heat 
exchanger by a centrifugal compressor and associated ducting. Steam 
is passed through the other side of the exchanger by the steam system. 
The air and steam systems are shown diagramatically in Fig. 4 and Fig. 5. 


The steam system is designed to give close control of the steam state 





within the core. Pressure and temperature instrumentation are pro- 
vided for measuring the air properties upstream and downstream of the 
core, the mass rate of air flow, the condition of the steam entering 
and leaving the core and the mass rate of excess steam and condensate 
leaving the core. 

The test apparatus is shown pictorially in Figs. l, 2, and 3. 

Air System. The air system is made up of 16 gage galvanized 
steel ducts joined with 1/4 inch steel flanges. The ducting prior to 
the test section has a maximum dimension of 12 inches square and the 
ducting after the test section is 14 inches in diameter. The lengths 
of ducting were designed in accordance with the ASME Power Test Code 
{7] for establishment of uniform flow and are shown in Fig. 4. 

The entrance to the test duct is on the exterior of the building 
to minimize thermal gradients in the air stream. The entrance is three 
feet square and is covered by a fine mesh screen to eliminate foreign 
matter from the air stream. The transition from three feet square to 
12 inches square is made in the entrance section, each side of which has 
the curvature of a quarter ellipse. After a five foot length of duct, 
a second transition is made from 12 inches square to the test core 
dimensions. This is followed by the forward instrument section which 
contains four thermocouple taps, a piezometer ring and two pitot tube 
taps for conducting vertical and horizontal velocity surveys. Follow- 
ing the test core is the after instrument section which contains nine 
thermocouple taps and a piezometer ring. A transition is then made to 
a 14 inch outside diameter duct at the end of which is a standard ASME 


square edged orifice. Piezometer rings are located one diameter 


10 





upstream of the orifice and one-half diameter downstream of the orifice. 
The duct section downstream of the orifice also contains a thermocouple 
tap. 

Air flow is regulated by a manually oparated double sliding plate 
valve located at the turbo-compressor inlet. The valve plates slide on 
Teflon runners and the drive screw is ball bearing mounted to eliminate 
maintenance problems. The double plates constrict the flow symmetrically 
and facilitate fine control of the air flow. Coarse control of the air 
flow is accomplished with the blast gate (discharge valve) of the cen- 
trifugal compressor. 

Air is induced through the test section ducting by a centrifugal 
compressor rated 6000 cfm at 2.5 psi back pressure. The compressor is 
driven by a three phase, 220 volt, 100 hp motor. The compressor dis- 
charges air through a 20 inch square duct to the exterior of the 
building. 

The air flow is metered by means of a standard ASME square edged 
orifice designed to the specifications given in the ASME Power Test 
Code (7). Four orifice plates made of 1/4 inch type 304 stainless steel 
were designed with the following orifice diameters. The preceeding duct 


has an inside diameter, D, of 13.875 inches. 


10.406 0.75 
6.244 (nce 
3.468 0225 
2.081 O.15 


These diameters were selected so that flow rates of from 250 lbs per hour 


to 22,000 lbs per hour can be measured with the ranges overlapping. The 


orifice plates are fitted with a centering grove on the bottom and a 


it 





scribe mark on the top so that they may be accurately positioned. 

The ducting downstream of the test section is insulated with two 
inch fiber glass insulation covered by aluminum foil, and the air duct 
is insulated from the test core by 1/8 inch Teflon gaskets. 

Steam System. Heating of the test core is accomplished by passing 
saturated steam through one side of the core. Before entering the core 
the steam is reduced to a low pressure, about six pSig, and a slightly 
superheated state, about five degrees superheat. In a superheated 
condition, the state of the steam can accurately be established by 
pressure and temperature. The few degrees of superheat are quickly 
lost in the test core and have negligible effect on the calculations. 

A considerable excess of steam is passed through the core to prevent 
a thick film of condensate from building up on the heat transfer sur- 
faces. 

Steam is supplied from the central steam heating system at 125 
psig in a saturated state. It is strained and then reduced to about 
30 psig through a 1 - 1/4 inch air operated, pilot controlled pressure 
reducer. At this point water is injected, if necessary, to desuperheat 
the steam and then it flows through a 1 - 1/4 inch centrifugal separator 
which removes all entrained moisture. The steam then passes through a 
second 1 - 1/4 inch air operated, pilot controlled pressure reducer 
which reduces it to about six psig. This second reduction slightly 
superheats the steam. Just prior to entry into the core the steam 
passes through a short transition section. To provide uniformity of 
flow in the test core, the transition section is filled with aluminum 


shavings held in place by a fine mesh stainless steel screen. Steam and 


2 





condensate from the bottom of the test core flow through a second 
transition section to a 1 = 1/4 inch centrifugal separator where the 
condensate is removed. A float type trap receives the condensate and 
discharges it at a realitively constant flow rate. The condensate is 
then subcooled in a small tap water counter flow heat exchanger and 
collected in a bucket for weighing. The exhaust steam flows from the 
Separator through a standard ASME square edged orifice with flange 
pressure taps and is discharged to the atmosphere. 

Two orifice plates were selected to measure the steam flow rates 
and have the following orifice diameters. The inside pipe diameter, 


Ds, is equal to 1.25 inches. 


d., 
ah (Gey) Pee De 
0.700 ~560 
0.890 ee 


The metering orifice was installed as specified in the ASME Power 
Test Code 9/4. 

Pressure and thermocouple taps were installed in both the upper 
and lower transition sections and a thermocouple tap was installed 
upstream of the steam orifice. The entire steam system is well insu- 
lated to minimize heat loss. 

Compressed air for actuation of the two steam pressure reducers 
is supplied by two air pressure regulators located on the instrument 
panel. This combination is designed to hold the downstream steam 
pressure to me 2 psig. An air line from the air supply of the first 


pressure reducer was fitted into the steam strainer clean-out plug so 


that the entire steam system can be blown dry upon completion of a 


iS 





series of test runs. This should minimize rusting when the system stands 


idle for long periods of time. 
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Pressure Instrumentation. Pressure instrumentation is provided 
for measuring , the gage pressure in the air duct upstream of the test 
core and upstream of the air orifice, the air pressure differential 
across the test core and air orifice, the steam header pressures and 
the pressure differential across the steam orifice. Air pressures are 
measured with well-type single leg water manometers and inclined draft 
gages. Steam pressures are measured with well type single leg mercury 
manometers. 

A one inch and three inch inclined draft gage and a 60 inch water 
manometer are used for measuring the air pressure differential across 
the test core. A 30 inch water manometer is used for measuring gage 
pressure upstream of the test core. The pressure tap downstream of the 
test core is located so that full pressure recovery of the air stream 
is attained. The air orifice pressure differential is measured with a 
three inch inclined draft gage and a 30 inch water manometer. A 60 inch 
water manometer is used to measure the gage pressure upstream of the 
air orifice. Each pressure tap on the air duct consists of four 1/16 
inch holes drilled symmetrically around the duct connected together by 
a piezometer ring of 1/4 inch copper tubing. One quarter inch plastic 
tubing is used to connect the piezometer rings to the manometers and 
draft gages and each manometer and draft gage is isolated by a brass 
toggle valve. 

Three 30 inch mercury manometers are used to measure the gage 
pressure of the top and bottom steam headers and the pressure differential 
across the steam orifice. Water is used as the intermediate fluid 


between the steam and mercury and a constant head of water is maintained 


20 





by using water pots. The change of level in the water pots is the same 
order of magnitude as the change in level in the mercury wells and these 
changes have insignificant effect on the pressure readings. The water 
pots are stainless steel cylinders three inches in diameter by three 
inches high and are located above the top of the manometers so that air 
pockets will not form in the water lines connecting them to the manometers. 
The connecting water lines are of 1/4 inch stainless steel tubing and 

each manometer is isolated by a stainless steel toggle valve. The water 
pots are connected to the steam headers and steam orifice by 5/8 inch 
copper tubing. The large diameter tubing was selected so that conden- 
sate can flow freely back into the steam system and not form pockets 

which would change the static head on the manometers. 

Temperature Instrumentation. All temperatures are measured with 
copper-constantan thermocouples. The thermocouples are in assemblies 
manufactured by Honeywell under the trade name MegapaK. The Megapak 
consists of a "head" for connecting the extension wires, a "Sheath" of 
1/8 inch stainless steel tubing with the insulated leads extending down 
through the center and the measuring junction at the end of the sheath. 
The thermocouples used for measuring the air stream temperatures are 
exposed,’ meaning that the junction is extended one sheath diameter 
past the end of the sheath. Those used for measuring the steam tempera- 
tures are "remote,'' meaning that the junction is one sheath diameter 
from the end of the sheath and the end of the sheat is sealed for pressure 
tightness. All extension wire is polyvinyl covered 24 gage copper-con- 
stantan. 


The thermocouple voltages are measured with a Honeywell 24 point 
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Multipoint recorder. The reference junction is an ice junction and is 
included in the measuring circuit. The recorder is a D.C. millivolt 
recorder with a scale of zero to ten millivolts. The smallest division 
on the chart paper is 1/10 of a millivolt and temperatures may be 
estimated to 1/100 of a millivalt, The recorder has several modes of 
operation. On the Multipoint setting it will balance and print each of 
the 24 points consecutively at a rate of about two seconds per point. 
On the Selecto=Print setting the recorder will balance and print at 
about two second intervals any number of the 24 points. On the Hold-on 
setting it continuously balances on any one of the 24 points but will 
not print. This setting may be used for obtaining a pen trace by 
attaching a clip-on pen to the print mechanism. The recorder’ should 
be adjusted for zero and span before use and the adjustment procedure 
is given in the operating procedures. 

The temperatures of the air upstream and downstream of the core, 
and the downstream of the orifice are recorded. Steam temperature in 
the top and bottom headers and upstream of the steam orifice are also 
recorded. To obtain an average temperature of air upstream and down- 
stream of the core a system of thermocouples is used. Upstream of the 
core this consists of three thermocouples placed on the horizontal center 
line. Downstream of the core nine thermocouples are arranged in a three 
by three grid with equal spacing between each of the thermocouples and 
the thermocouples and the wall. In this manner small variations in 
temperature across the stream can be averaged out. All other tempera- 
tures are measured by a single thermocouple. Each penetration through 


the duct for a thermocouple is fitted with a Swagelok compression fitting 
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so that the penetration length can be adjusted. 
The thermocouples are numbered as shown in Figs. 4 and 5 and the 
Same number is printed beside the temperature reading on the multipoint 


recorder chart. 


4, Experimental Method and Equations 

General. The experimental method and equations used are essentially 
those outlined by Kays [3] for a similar facility. These equations are 
suited for use with tables and graphs. The equations outlined herein are 
suited for use in a digital computer program which will accept raw test 
data and produce finished results. The equations used in preparation of 
original graphs and tables were used if available and approximate 
equations were written for those not available. Graphs of the approxi- 
mations are shown in the appendix and test data should be checked for 
applicability of the approximations. 


Mass Rate of Air Flow. The mass rate of air flow, mn. through an 





ASME square edged orifice. is given in the Power Test Code [7] as: 


m = 359 CFd°F,YVAR,® Ib, “nr 


where C 


coefficient of discharge 


F = velocity of approach factor 

d = orifice diameter, in. 

Le = orifice thermal expansion factor 

x = net expansion factor 

AP, = pressure differential across the orifice, in. H50 

Ro = specific weight of the air upstream of the orifice, 
1b,,/£t3 
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The ratio of the orifice diameter to the duct diameter is: 


= <@d 
a 
where D = duct diameter 


13.875 tn. 


and the velocity of approach factor is: 


The absolute duct pressure upstream of the orifice is given by: 


Py = Pi (0-490) — Py (0.036) psia 


where P, = barometric pressure, in. Hg 
Eo — welllcE (pPlessure Upsieream or fhe Oririce, sin 150 
(duct pressure is below atmospheric DreeeneeN 


A dimensionless ratio of orifice pressure drop to orifice pressure 


is: _ AP, (0.036)) 
x = PR, 


and the net expansion factor is then: 
Y = 1.0- (S)coa + 0.35%) 


The approximate equation for the orifice thermal expansion factor 


for type 304 stainless steel gives: 


eS =e lea + (to a 68)(1.85 x iO”) 


The coefficient of discharge is an unknown in the equation but may 
be quickly iterated using an equation suggested by Murdock [6]. The 


equation is dependant upon the Reynolds number in the orifice, Nred>? and is: 


4\\~% 
C= Cy + ac fio" 


N Red 
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The coefficients @, and AC for applicable 8's from Murdock [6] and a 


suggested first iteration C are: 


B Ge AC C 
ni 0.59446 0.00945 0.5975 
25 0.59483 0.01037 0.5966 
“AS 0.59863 0.01543 0.6014 
is 0.60480 0.05448 0.6128 


The density of the air upstream of the orifice may he found from 
the equation of state of a perfect gas modified by a humidity correction 


factor which gives: 


_ 1442 Xm 
Go= 5334 7. 


53.34 To 
where Po = duce pressure Wpeeream Of the oritrice, psia 
Too se LenPersbt meso —ealesuns Eeeamo. athe somirtiec. OR 
X, ™= humidity correction factor for density as given by 


Kays and London [4] 


(+H 
(1 + 1.607 H) 
and H = humidity ratio, lbs, Hj0/1b, dry air 


The dynamic viscosity, AA,, for air at one atmosphere pressure is 


3/2 
Wo = ©, OOMib2 let 
Tok + 110.4 


given in reference (1) as: 


where Tok = temperature of air in duct, 7K 


The orifice Reynolds number may be found from: 


N _ 15.28 mh 
Red = Mod 
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Heat Transfer Calculations. The unit conductance for thermal con- 





vection heat transfer as given by equation one in the theory section is: 


| fi 
h =— m_ 1 aa a +- ee WY Ons 
3 Sa 3 
No} Uy "S| Ne Ag he 2A kw 
where he = unit conductance for thermal convection heat 
Clanceereromrnowaim.) 6 G0) (in ft? OF) 
Yo = total surface effectiveness for the air side 
U = overall unit conductance, BTU/(hr ft2 °F) 
ne = total heat transfer area on the air side, £t2 
or = total surface effectiveness for the steam side 
i = total heat transfer area. on the steam side, ft? 
hg = unit conduction for thermal convection on the 
steam side, BTU/(hr ft2 °F) 
Q., = thickness of the wall separting the steam side 
from the air side, in. 
Awa = area of the wall on the air side, ft? 


thermal conductivity of the wall, BIU/(hr ft? OF/ft) 


& 


The film coefficient for convection on the steam side, Ne is unknown 
and must be estimated. Any reasonable estimate is considered adequate 
since a 50 percent error in the estimate results in about two percent 


error in hg. A value of hg equal to 2000 BTU/(hr ft2 °F) has been used. 


The specific heat of dry air may be considered constant and is assumed 


pombe Coqa = 0-241 BTU (Ibm °F) 


The humidity correction factor for the specific heat as given by Kays 


and London [4] is: 


, [= IS) Sr. 
c* 1 +H 
where H = humidity ratio, lb , H50/1b, dry air 
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and the corrected specific heat is then: 


The bulk average temperature in the test core is: 


tye o 
tp 


where Le) air temperature upstream of thetest core, OF 


t9 air temperature downstream of the test core, °F 


and the film average temperature is then: 


tptt, ° 
tp= “bos *F 


where ts = the saturation temperature corresponding to the 
arithmetic average steam pressure in the core, °F 


The number of heat Transfer Units may be found from the temperature 
NTU relation as given in the theory section, yielding: 


ts- 


NTU = &n 
Us a be 


and the overall unit conductance is then: 
& 
eam ee o 
_ CpmN’* atu tt® °F) 
Sea 
a 
For most plate-fin type geometries the equation for the temperature 
effectiveness of a straight fin with constant conduction cross section 


may be used to a good approximation. [3]. The fin temperature effective- 


ness for the steam side, Ngg» is: 


tanh (ms 257 12) 
(m9, 12) 


where £, = effective fin length for the steam side and is 
equal to one-half the total fin height, in. 


Neg 


24 Ns —| 
Me = ciara fy 
? Kes Ss 


Ze 





and 


2000 BTU/(hr £t2 °F) as previously estimated. 


thermal conductance for the fin on the steam 
side, BTU/(hr ft? °F/ft). 


fin thickness on the steam side, in. 


The overall temperature effectiveness for the steam side, Ng, is then: 


where 


The fin temperature 


where 


and 


The overall temperature effectiveness for the air side, n 


Awe + Nea Afa 


where 


Ns 
As 


Ags 


As 


La 


— 
—_— 


Awe ate Nee Ats 
Ag 


wall area on the steam side, ft 


fin area on the steam side, ft? 


total heat transfer area on the steam side, ft 


effectiveness for the air side, ) fa; es’ 


Ge 


tanh (m5 t, /2) 
Mala 12 


effective fin length for the air side and is equal 
to one-half the total fin height, in. 


24 ha ad 
e Jot ft 
Me Key tg 


No = 


A 
wa 


Ata 


Aa 


thermal conductance for’the fin on the air side, 
BTU/(hr ft? OF/ft). 


fin thickness on the air side, in. 
film coefficient for convection on the air side. 
It may be taken equal to U for the first iteration 


BTU/(hr ft2 OF). 


o» is given by: 


Ag 
wall area on the air side, ft? 
fin area on the air side, ft 


total heat transfer area on the air side, £t2 
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The air mass velocity in the core, G, is given by: 





° 2 
. 2h ip. faunr ft 
or Ac of ' 


where A 


a the free flow area of the core on the air side, ft2 


Stanton number in the core can now be evaluated as: 


es 
"St Sey 


The dynamic viscosity evaluated at the film average temperature and 


one atmosphere pressure as given by Hilsenrath [1] is: 


3/2 
Ue = 0.003527 Tek 1D /(he ful 
T Te. + HO.4 


fk = film average temperature, °K 


where T 


The dynamic viscosity evaluated at the bulk average temperature is: 


3/2. 
0.003527 Tbk Ib Reuse 
= tc 
Mom T+ 11004 m/ (re £¢) 
where Ti, = bulk average temperature, OK 


Ther thermal conductivity of air at the film average temperature 


and one atmosphere pressure as set forth by Hilsenrath [1] is: 


V2 
i = 0.001529 Tk BTU 
TT, 4 245.4 i Oh eae oe 
Te X10°7 Tk 
Where ere = film average temperature, °K 


The Reynolds numbers evaluated at film average temperature and bulk 


average temperature are evaluated by: 


4rhG 


N Ref 7 Mr 
and 
4r4G 
2) eee 
NReb ar Mb 
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where hydraulic radius r, is defined by Kays and London [3] to be based 
on the free flow area and takes the form 
Nell 
rh — 
Ag 


where L 





Ft 


core flow length, ft 


Prantdl number in the core is: 


- Ais 
= 


where the fluid properties are evaluated at the film average temperature 
in the test core. 


The Colburn-j heat transfer modulus is then: 


Se 2/3 
J = No Np- 


which concludes the heat transfer calculations. 
Friction Factor Calculations. The equation for the core friction 


factor as developed in equation two in the theory section may be written 


as: rh Om AP. (4.3255 x ia) Ke —(| AF g-¢) 
SS SSeS a> ~~ 
L Ge | 


fara 
Ke tl + 0-5 (+4 Dg } 
Qo 


where r, — ene hydraulic radius, ft 
Om = mean density in the core Ib,,/ ft? 
L = core flow length, ft 
APc = pressure drop across core, in. H,0 
G = core mass velocity, 1b,/hr fr2 
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Dq 


contraction coefficient 

the ratio of free flow area to frontal_area 
density of the air upstream of core, lb/ft? 
expansion coefficient 

8 


density of the air downstream of core, 1b,/ft 


friction factor for the duct downstream of the 
core 


length of duct from the core to the downstream 
pressure tap, in. 


hydraulic diameter of the duct downstream of 
core, in. 


The duct friction factor, fd, is assumed constant at .015, a value 


corresponding to a nearly smooth pipe for a Reyholds number range of 


100,000 to 500,000 as may be expected in the downstream duct. 


The mean density in the case of heating is given by: 


where 


P9 


TO, + Tol. 


3 
: 1s - (= 72) lbp, /Ft 
NTU 





upstream temperature, OR 
downstream temperature, OR 


Saturation temperature corresponding to the average 
steam pressure in the core °R 


upstream density 


0.097318 X,P 3 
ea eee 
T, mM awit 


upstream pressure, in. H,0 absolute 


humidity correction factor as previously given 


downstream density 


- 0.09318 Xm P2_ ib, Ar+3 


T2 


Sl 





B, = downstream pressure, in. H50 absolute 
In the isothermal case, the mean density is equal to the arithmetic 


average density and since the temperature remains constant 


_ 0.097318 X,,(P, + Pe) Ib, /r 43 


Om = Cave = 
Zily 
The equations for the contraction coefficient, Kc, and the expansion 


coefficient, Ke, are developed by Kays [2] and are given by 


_. 


2 
and Ce 
e le 
where i see LOctiyedi1 Str rpilmions coer ilcienir 
g = be / Afr. 
C. = jet contraction ratio 


The jet contraction ratio is approximated as 

C. = 0.611 +0.45 + 0.344 06>! 
Thesveloetey-distribucion coetficient «Ka is dependent upon both core 
geometry and film Reynolds, Np of’ The usual geometries used in plate 


fin type construction are: 


square at | oe 
triangular VPS 


For laminar flow, Npef less than 2000, Kays [2] gives the values 
square Ka Ne eee 
triangular Kg = Kq, = 1.43 


and for turbulent flow 





square Ka = 1 oy (Ka. - 1) 


i = + - ns 
triangular Kq Kq 1. 1229 (Ka. 1) 
where the velocity-distribution coefficient for circular tubes, Kde > is 
given by: 


Kd¢ = 109068 (4fm) +0:0584 -+4F, +1" 


and —O.2 
f.1 = ©0429 Noor 


If the boundary layer is continually tripped, as in the case of louvered 
tin or off-set rin construction, film Reynolds number should be taken 
equal to infinity which ygelds Kd = 1.0 in all cases. 
Heat Balance Calculations. The heat transfer rate to the air, 
Qwair, is: 
Q air = m Cp Ge5e— t1) BIU/hr 
The total heat given up by the steam, Q steam, is equal to the heat 
given = by the excess steam as it passes through the core plus the heat 
given up by the condensate in the core. 
Q steam = § (hy ~ hog) + % (hy - hof) BTU/hr 
where Ss = Mass rate of flow of excess steam, lbp/hr 


mass rate of flow of the condensate, 1b,/hr 


= 
Q 
Il 


hj = enthalpy of the steam at the top header pressure and 
temperature BTU/1b_, 


Doe = 


2 enthalpy of saturated steam at the bottom header 


pressure, BIU/1b,, 


ho¢= enthalpy of saturated liquid at the bottom header 
pressure, BTU/1b_, 


o 
The mass rate of flow of the excess steam, Ss, as measured by a square 


edged, standard ASME orifice with flange pressure taps is found from the 
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applicable equation given in the ASME Power Test Code [7]: 


§ = 359 CyFodoFacvs,/ Pe Ib, “hr 


Vg 
where ce = coefficient of discharge 
He = velocity of approach factor 
d. = orifice diameter, in.. 
FAs = Orifice thermal expansion factor 
Mie = net expansion factor | 
Ap. = pressure differential across the steam orifice, 
Tis Ho0 
Vs = specific volume of the steam passing through the 


orifice, fc aes 


The velocity of approach factor, Fg,, is: 


Bl le 
s Vi- 6 


iT 


where 


p,= ae 
S De 
and Ds = steam pipe diameter 


1.25 ime 
The orifice thermal expansion factor is taken as constant corres- 


ponding to a steam temperature of 220 OF. For type 304 stainless steel 


a = 1.003 

S 

The net expansion factor for steam, Yg, is: 

Ys =1— (0.41 + 0.3563) (%s/103) 
xg= APs/Ps 


The upstream steam pressure, P 


where 


g» is taken as the saturation pressure 


corresponding to the upstream temperature and is given by the approximate 
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equation 


Ps = 16.5 + (to —218)(0.376) PSta 


The specific volume is given by the approximate equation: 


a 3 
Vs = 24-01—0.4(tso —218) ft/ Ib 
The coefficient of discharge, Cos is iterated by the same method 
that was used for the air orifice. The values of ee and AC for applicable 


values of Bg, are: 


Be die Cy = Cos AC = AC. 
56 . /00 in. ~6085 -025 
Pel .890 in. . 6090 ~041 
and then 


\/2 
a (08. 
esd 


The steam orifice Reynolds number, Npagq, is: 


' _ 0.004244 § 
Resd Ms ds 


and the dynamic viscosity is given by the approximate equation: 


= | =8\ _!f 
Ms = 8925 x 10 © + (tg9— 200\(1-45 X JO ) ft pec 


The enthalpy of the superheated steam in the top header is approxi- 
mated by two equations, the first giving the variation with tempgrdture 
at a constant pressure, P.4 = 18.0 psia and the second equation the 


variation with pressure. These equations are: 


/ 
hey = IT50s0 4 Cte 2300) (15) BIU/1b,, 
where Bee = enthalpy at 18.0 psia: 
tsl = top header steam temperatrue, °F 


a2 





thence hey ha = CE ae tomes s5) BLU/ lb, 


where Pol steam pressure in top header, psia 
The enthalpy of the saturated steam in the bottom header is found 


from the approximate equation: 


hg2 1155.3 + (Poy - 19.0) (.934) BIU/1by 


where P32 stream pressure in bottom header, psia 
and the enthalpy of the condensate in the bottom header is: 
ie piel. 4 tee 19 O72 64) BTU/1b,, 


The mass rate of flow of the condensate, w is found by timing 


ce? 


the condensate collected. 


5. Evaluation 

A steady state steam to air test was conducted on a compact heat 
exchanger manufactured by Hamilton Standard Corporation. The test core 
had the approximate dimensions of a six inch cube. The construction was 
the plane plate-fin type with flow passages between fins being approxi- 
mately square. One side of the core was "straight,'' meaning that the flow 
passage was straight for the length of the exchanger. The other side was 
"ruffled," meaning that the flow passage was approximately a cosine curve. 
The plate spacing was 0.156 inches and the fin spacing was 18 fins per 
inch for both sides. All plate and fin material was 0.004 inch stainless 
steel. A set of "standard" curves for Colburn-j and friction factor vs. 
Reynolds number were supplied with the test core. The "standard" curves 
were for plate and fin spacing similar to the core tested. 

The test was performed using the straight side of the core for air 
and the ruffled side for steam. The various compgnents of the facility, 
the air system, steam system and instrumentation were checked out and 


proved to be very satisfactory. There was a problem in defining flow and 
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heat transfer areas for the exchanger because structural stiffeners had 
been placed at the end of each plate-fin passage. The stiffeners were 
placed at an angle across the flow path and thereby restricted the 
flow in the peripheral flow passages. This restriction would tend to 
increase with increasing flow rates. The test results based upon the 
minimum free flow area, which neglects the peripheral flow passages 
affected by the stiffeners, was plotted on the standard curves as 
shown in Fig. 6. 

The next test should be on a heat exchanger core with previously 


determined characteristics and with given or definable geometry. 
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Figure 6, Colburn-j and friction factor 
vs. test core Reynolds number. 
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APPENDIX I 
OPERATING PROCEDURE 

General. The preliminary check of instrumentation for both the 
hot core and cold core tests is the same. 

The multipoint recorder should be allowed to warm up two hours 
prior to calibration. To calibrate the instrument, set the print 
mechanism to “Hold On" indicating any one of the points not used for a 
thermocouple. This open circuits the thermocouple measuring circuit. 
Place a D.C. millivolt source across the terminals at the back of the 
recorder and set the source to one millivolt. Adjust the recorder zero 
so that one millivolt is indicated. Set the source to nine millivolts 
and adjust the span so that nine millivolts are indicated. Repeat the 
zero and span adjustments until the recorder is calibrated. 

All manometers should be adjusted for zero with the isolation valves 
open prior to lighting off the air or steam systems. The two steam 
header mercury manometers must be zeroed at about three inches due to the 
head of water between the manometers and the water pots. The water lines 
between the mercury manometers and the water pots are susceptible to the 


formation of air pockets which will cause erroneous readings. 
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Cold Core Test. In order to assure an even distribution of data 
points, the approximate pressure differential across the air orifice 
should be calculated for each run. Calculate the mass rate of air 
flow at a core Reynolds number of 500 and 10,000 and plot these two 
points on a log log plot. Divide the straight line connecting the points 
into equal increments, the number depending upon the desired number of 
runs. The orifice pressure differential may then be read from Fig. /7 
for each mass rate of air flow. 

Open all manometer valves and adjust for zero. Close all manometer 
valves prior to lighting off. Set the multipoint recorder to print the 
thermocouples in use by raising the corresponding buttons on "Selecto- 
Print.'' Turn the chart motor on and select the desired speed. Close 
the compressor blast gate (discharge valve) and energize the compressor 
motor. Open the blast gate to desired setting. Adjust the sliding 
plate valve for the desired air orifice pressure differential and observe 
instrumentation until steady state obtained. 

Record the following data for each run: 

1. Atmospheric pressure (0.001 in, Hg) 

2. Atmospheric dry bulb temperature (0.1 °F) 

3. Atmospheric wet bulb temperature (0.1 mE) 

4. Air orifice size 

5. Orifice upstream pressure (0.01 in, H,0) 


6. Orifice pressure differential (0.001 in, H,0 for APo <3", and 
0.01 in.H50 for APS >3") 


7. Core upstream pressure (0.01 in, H,0) 


8. Core pressure differential (0.001 in,H,0 for APc <3", and 0.01 
in,H,0 for APc >3"). 
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9. Orifice air temperature (0.01 millivolt) 

10. Core upstream air temperature (0.01 millivolt) 

Items five through ten are recorded three times in the same sequence 
and averaged on the data sheet. 

Hot Core Test. The instrumentation checks and the procedures for 
lighting off the air system are the same as for the Cold Core Test. 

The steam system is activated by the following procedures. Obtain 
100 psig low pressure air to the reducer central valves. Turn on the 
tap water supply for the condensate cooler and the desuperheater. Crack 
the desuperheater supply valve. Open the by-pass valves for the steam 
reducers and crack the steam stop valve. Allow the system to warm up 
slowly. When warm-up is completed close the by-pass valves and open the 
steam stop. Charge the first steam reducer to about 30 psig and the 
second to about six psig. Obtain a steam state in the top steam header 
of about six psig with five to ten degrees superheat by adjusting the 
first steam reducer and the desuperheater. The excess steam flow rate 
should be at least five times the condensate flow rate. Two steam 
orifice plates are provided; the larger to be used for high air flow rates. 
Condensate should be collected for at least five minutes at high air flow 
rates and for at least ten minutes at low air flow rates. 

Steady state must be achieved before each run. This will require 
about one hour prior to the first run and 15 minutes between runs. Steady 
downstream core temperatures are the best indication of steady state. 

The following data is recorded for each run: 

1. Atmospheric pressure (0.001 in. Hg) 


2. Atmospheric dry bulb temperature (0.1 °F) 
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10. 


rs 


2% 


i 


14. 


i. 


16. 


1/7. 


18. 


19. 


20. 


Atmospheric wet bulb temperature (0.1 °F) 

Air orifice size 

Steam orifice size 

Top steam header pressure (0.01 in. Hg) 

Bottom steam header pressure (0.01 in. Hg) 

Pressure differential across the steam orifice (0.01 in. Hg 
under H90) 


Core upstream pressure (0.01 in. H50) 


Core pressure differential (0.001 in. Hj0 for APc< 3", and 0.01 
in, H 0 for APc> 3") 


Orifice pressure differential (0.001 in. H,0 for APo< 3", and 
0.01 in.H»90 for APo >3") 


Orifice upstream pressure (0.01 in. H,0) 

Core upstream air temperature (0.01 millivolt) 
Core downstream air temperature (0.01 millivolt) 
Orifice air temperature (0.01 millivolt) 

Top steam header temperature (0.01 millivolt) 
Bottom steam header temperature (0.01 millivolt) 
Steam orifice temperature (0.01 millivolt) 
Condensate weight (0.1 1b) 


Condensate collection time (l. sec) 


Record items 6 through 18 three times in the same sequence and 


average on the data sheet. 


The humidity ratio, H, may be found from a standard psychometric 


chart and convert temperatures to OF (OL). 
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Figure 7. Mass rate of air flow vs. 
orifice pressure differential. 
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APPENDIX II 
DIGITAL COMPUTER PROGRAM FOR DATA REDUCTION 

The program is designed to reduce data from both the hof core test 
and the cold core test. The program written in FORTRAN 60 and is shown 
on pages 42through 9/ and the output for the hot core test is shown on 
page 58 and that of the cold core test on page 59 . The variables used 
in the program are defined in the program glossary. Variables not appear- 
ing in the glossary are internal to the program and are used to change 
dimension or to define special groupings to simplify programing. 

The index used to reduce data from the hot core test or the cold 


core test is given by: 


Program Index 


IHORC 
Hot core test 1 
Cold core test iz 


The program is indexed to accept the various fin geometries normally 


used in plate fin construction and the indices: are as follows: 


. Fin Geometry Index 


‘ 


ISQTRIN 
Square fin 1 
Triangular fin 2 
Louvered or off-set fin 3 


(Ne = infinity) 
The indices used for the fluid metering orifices are: 
AipeOririce Air orifice Index 


BETA d (in.) IB 
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215 10.406 4 


45 6.244 3 
oe) 3.468 Z 
5 2.081 l 
Steam Orifice Steam Orifice Index 
BS dsms(ini. } IBS 
/1 0.890 Z 
0 0.700 1 


The indices are all entered as fixed point constants. The number 
of runs, NOR, and the run number, NR, are entered as fixed point con- 
stants. All fixed point constants must be right justified in their 
data fields. 

The constants for a given core such as areas, lengths, etc. are 
entered as input data only once and are the first three data cards 
beginning with AC and ending with DD. The fixed point constants for a 
particular series of runs, NOR, IHORC, ISQTRIN are entered on the fourth 
data card. The fifth through the seventh data cards are for run one, 
eighth through tenth for run two and so on. The first data card for a 
run contains the run fixed point constants, NR, IB, IBS. The second and 
third data cards for a run contain the run data beginning with H and 
ending with TC. When entering data for the cold core test, only the 
applicable data must be entered. All other data fields may be left blank. 

The dimensions of input data must be those given in the program 


glossary and are the same as the raw data. 
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Program 
Symbol 


AA 
AC 
ALFA 
AFA 
AFR 
AFS 
AKFM 
AS 
AWA 
AWS 
BETA 


BS 


CC 
CO 
COLJ 
COS 
Cr 
CRKD 
CS 
DC 
DCS 


DD 


PROGRAM GLOSSARY 


Program 
Unit 


et 
ft 


ft? 


£t- 


ft? 
BIU/(hr f£t2 °F/ft) 
ft? 
ft? 


ft? 


Blu bE) 


oe 
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Nomenclature 
Symbol 


Aa 


Ag 





DL 


DPC 


DPO 


DPS 


DSO 


ERROR 


FA 


FAS 


FB 


FBS 


FC 


FD 


FLA 


FLL 


FLS 


FKA 


FKS 


FM 


HG 


HS 


HS1 


HS18 


Hoz 


dni. 

in. H50 

in. H90 

in. Hg under water 


iti. 


il. 
et 
Enis 
BTU/(hr f£t20F/f£t) 


BRU Chir ote2eOr/ tt) 


lbs/(hr ft7) 
1b/1b 

BTU/ (hr f£t2°F) 
BTU/1b__ 
BTU/(hr £t2°F) 
BTU/1b_ 
BIU/1b,, 


BTU/1b,, 
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IB 
IBS 
ITHORC 


ISQTRIN 


NOR 
NR 
OD 
Pl 
P2 
PB 
PO 
PR23 
PS1 
PS2 
PSO 
Q AIR 
QSTM 
REB 
REFM 
REO 
RESO 
RH 
RO1 
RO2 


ROAVE 


Index for air orifice diameter 
Index for steam orifice diameter 
Index for Hot or Cold core test 


Index for fin geometry, square, triangular 
or infinite Reynolds number 


Number of runs 


Run number 


in. d 

in. H90 Ps 

in. H90 Po 

in. Hg Er 

in. H90 le 
No.2! 

ine He Pol 

in. Hg Po2 

psia Eas 

BTU/hr Q@aar 

BTU/hr Qsteam 
NReb 
NRef 
NRed 
NResd 

ie ae 

Lb! ,p/ft> 01 

Lb. 4,/£t> 02 

lb m/ft? oe 
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ROM 
ROO 
SDOT 
ST 
SQKD 
Tl 
Te 
TB 


TC 


TFA 
TFM 
TFS 
TO 
TS 
ESL 
ao 
TSO 
TRKD 


TW 


UB 


UC 


UFM 
UO 


USO 


| 3 
eee 


Pm 
Lb: ,/ft? pn 
lb m/ hr S 
Nst 
Ka 
S 
or ty 
= t5 
oe th 
min time of condensate 
collection 
ite Sa 
oe fer 
in Ss 
5 GA 
or ts 
“18 Cs1 
=e ay, 
OF Cso 
Ka, 
in’. X 
BTU/(hr £¢°F) U 
1b: {,/ Chr ft) Mb 
lovecat. £e) “in core for cold core 
Cest 
lb a Ghiew it) ak 
Lio ie) 0 
lb m/ (sec Et) LPS 
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vs ft3/1bm 


WC Lb. 


WCH iba / he 
WDOT Lb‘ p/hr 


WK BTU/ (hr £t2°F/£t) 


XC 


SKD 


XNTU 


XS 


YKE 
YNFA 
YNFS 


YNO 


¥S 


ZKC 


Dil 


Vs 


Weight of-condensate 
collected 


o 
We 


6 
m 
ky 
x 


X 
Cc 


Kq 


Xm 


Mg 


m 
S 


NTU 


Xs 





i 
2 


3 
4 


5 
6 


7 
8 


9 


20 


ag 


22 


28 


25 


PROGRAM SSHEAT 

OREAD 1,AC,AFR,AA,AWA,AFA,AS,AWSsAFS,FLL> TFS, TFA,TW, 
IFLA,FLS,»FKA»FKS,WK,FD.DL, DD 

FORMAT (8F 10.0) 

READ 2,NOR,IHORC,ISOTRIN 

FORMAT( 3110) 

GO TO (3.5). ITHORC 


PRINT 4 

OFORMAT(CLHL 25X%,32H HEAT TRANSFER AND FRICTION DATA,/S// 

1106H RUN MDOT G HUMID T1 Pi T2 pPc=OrsésT$ 

2H BTU/ COLU F FACT NRF NR B  ERROR,// 

S107H LB/HR LB/HR FT2 LB/LB F PSIA F PSI FO 
4R FT2 F PERCENT»//) 

60 wor7 

PRINT 6 

OFORMAT(1H1 25X,25H ISOTHERMAL FRICTION DATA,/// 

160H RUN MDOT G HUMID T1 P1 DPC oF FACT NR,/S/ 
246H LB/HR LB/HR FT2 LB/LB F PSIA PS]s//) 

READ 8,NR,18,1BS 


FORMAT (3110) 
ie SmdyOD,DSO,BETA,8S,TO,11,12,1S4,1S¢,1S50, 
1PB,PS1,PS2,DPS,P1,DPC,P0,DP0,WC,TC 
FORMAT (8F10.0) 
PS1A=.490*(PB+PS1) 
PS2A=.490%(PB+PS2) 
PSAVE=(PSLA+PS2A) 7/2. 
TS=220.+(PSAVE~17.1)%2,.66 
P2=P1+DPC 
TOR=T0+459.7 
TOK=TOR*5./9. 
UO=.,005527*TOKee1,5/(TOK+110,4) 
GO TO (20,21,22,235 ),I1B ae 
CO=.59446 
pc=.00945 
C1=.5975 
OD=10- 406 
BE Tae. 75 
GO WO 25 
co=,.59488 
ie = OOS 7 
C1=,5966 
OD=6.244 
BETA=.45 
GO TO 25 
CO=.59868 
pDc=,01543 
C1=.6014 
OD=3,468 
BETA=.25 
gion 1025 
COo=.60480 
nc=.05448 
C1=.6128 
Oop=2.081 
BETA=.15 
C=C 
B4=BETA**4 
POA=PB*,490~-P0*,03561 
X=(DPO*,0361)/P0A 
Y=1.-66(6.4914+, 55484) #xX)/1.4) 
XM=(1+*H)/(1.+1,607*H) 
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26 


29 


50 


40 


ai84 
42 


43 
44 


45 


ROO=PO0A*144.*#XM/(53.35*TOR) 
Me (10-66.)12,12,11 
FA=1.+(T0-68.)*1,85E-5 
G66 10 43 
PA 1 
FB=1./(1.-B4) **,9 
DO 26 I=1,3 ; 
WDOT=H=359.*C#¥FBeFAxY*ODeOD*eSQRTF (DPO*ROO) 
REOQ=15,28 *WDOT/(U0O*OD) 
C=CO+DC*SQRTF(10000./RE0) 
CONTINUE 
GO TO (29,40).IHORC 
TB=(T1+T2)/2. 
TFM=(TB+TS)/2. 
TBK=(18+459.7)*5./9. 
TFMK=(TEM+4959.7)4%9./9., 
UB=.003527*TBKweel .5/(TBK+110.49) 
UFM=.003527*TFMKee1.5/(TFMK+110.9) 
AKEM= OOL529"TFMK HR. S/ (1.4295. 4/ 0 TF MK10.e8(12./TF MK) )) 
RH=AC*FLL/AA 
G=WDOT/AC 
REF M=4.*RH*G/UFM 
REB=4.*RH*G/UB 
XC=(Le+1-915eH)/(1.,4H) 
BP= XC*,.241 
MN TPemOGr €(1S-14).7(1S-T2)) 
U=CP*xWDOT*XNTU/SAA 
HS=2000, 
XMS=SQRTF (29. *eHS/(FKS*TFS)) 
YNFS=(TANHCXMS*FLS/12-))/(XMS*FLS/12.) 
YNS=(AWS+YNFS*AFS)/AS 
HA = U 
DO $0 l=1,8 
XMA=SQRTE (24. *HA/S(FKA®TFA)) 
YNFAS(TANHCXMA®FLA/12-.) )/(XMARFLA/1I2.) 
YNO=CAWA+YNFA®AFA)/SAA 
WAEL A CYNO#((2./0) MAA# (1. / CYNS#AS#HS) +TW/(AWABWK RLS.) DD) 
PR23=(UFM*CP/AKFM) **.66/ 
ST=HA/(G*CP) 
COLUJ= ST*PR23 
ALFA= AC/AFR S 
P4A= PBei3.571- Pi s 
P2A= PBeiS,571- pe 
TLA= T1+459.7 
T2A= T2+459.7 
ROL=-097SLE*XMHPLASTIA 
ROZ=.0973518*XM*P2A/ST2A 
GO TO (41,42), 1HORC 
ROMS. S¥((TLARROL+T2ZA®ROZ) /(TS4459.7-((T2A*TIAY/XNTU))) 
GO 10 436 
ROAVE = .097318"XMe(PLA+P2ZA)/(2."TLA) 
ROM=ROAVE 
G=KDO017AC 
RH=AC¥FLL/AA 
TiK=(714+459-7)%*5./9. 
UC=.003527*11Kee1,.5/ (TL K+110,.4) 
REC=4.*RH*G/UC 
REFM=REC 
GO TO (50.45,44),ITSQTRIN 
XKD=1> 
GO TO 60 
IF (REFM"2000.)46,497,4/7 
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46 


47 


48 
50 
22 
53 
54 
20 
56 
60 


70 
75 


76 


77 


80 


TRKD=1.438 

GO TO 48 

FM =.0499*REFM**(-,2) 
CRKD=1.09068%4,%FM+.05884eSQORTF(4. "FMI +1, 
TRKD=1.+1.29%(CRKD-1.) 

XKD=TRKD 

GO TO 60 

IF CREF M -2000.)51,53,58 

SQKD = 1.389 

GO TO 56 

FM =.049*REFM**(-.2) 
CRKD=1.0906844,«F M+, 09884eSQORTF(4.*FM) +1, 
SQKD = 1.+1-17*(CRKD-1.) 

XKD=SQKD 

CC =.611+.045*ALFA +. 5449" (ALFA®*5.7) 
YKE=1-72.*XKD*ALFA+ALFA®H2, 
ZKC=(1.72-*CC+(CC#¥CC) *( (2. *XKD)-1.))/°(00*%CC) 
A=(ZKC-(1.+ALFA®ALFA))/ROL 
B=(YKE+1.+ALFAwALFA®(1.+(4.*FD*DL/DD)))/RO2 
FC= (RH*ROM/FLL) *((64.$255E+9*DPC/(G*G))-A~B) 
PTAA = PiA*.0561 

DPCA=DPC*,.0561 

GO TO (70,104),I1HORC 

Geo 10 (75,76),]BS 

Cos=.6090 

DCS=.0410 

DSO0O=.89 

BS=./71 

GO To 77 

COS=.6085 

DCS=.0250 

DSO=.70 

BS=.56 

CS=,61 

BS4=BS**4, 

FBS=1./(1.°BS4) 4.5 

FAS=1.008 

DPSO=.454«DPS - 
PSO=17.1+(TSO-220.*, 3576 

XS=DPSO/PSO 

¥S=4.70 2494+. 355"BS4)*(XS/1.5) 
DPSW=DPS0/.0361 

VS=24.01-(TS0-218.)*.4 
USO0=8.25E-6+(TS0-200.)*1.45E-8 

DO 80 15=1,3 
SDOT=359.*CSeFSBeFAS*eYS*DSO*DSO*SQRTF CDPSW/VS) 
RESO=.004244*SDOT/(USO*DSO) 
CS=COS+DCS*SQRTF(10000./RESO) 

CONTINUE 

HS18=1158.0+(17S1-250.)*.5 
HS1=HS18-(PS1i1-18.0)*.55 
HS2=1155.3+(PS2A-19.)%(.954) 
HC=195.4+(PS2A-19,)*2.64 

WCH=WC*60./7C 

QAIRaWDOT*CPx(T2-T1) 

QSTM=WCH* (HS1-HC) +SDOT*(HS1-HS2) 
ERROR=(QAIR-QSTM)*#100./0AIR 
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99 PRINT 400,NR,WDOT,Gs»H»T1,PLAA»T2,DPCA,TS,HA,COLJ,FC,REFM,REB, ERROR 
LOQOFORMATC IS, 2X oF 6.1, 2X0 F749 2X oF 564s 2X oF 4 ole 2XoF Se 2r2XsF Suir eXs Fh Gee, 


104 


12X pF 501) 2X pF 4ed, SX pF 605s 2X FO. dr eXeF 6.4, 2Xo Fb. Ls, 2XsF5.1s5//) 


GO TO 110 
PRINT 105,NR,WDOT,.G,H,T1,PLAA,DPCA,FC»REFM 
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LOSOFORMAT (135.2% 0F6.102X oF 7 odo 2X F504 2X oF 4dr QXoF 52 r2XsF 420 2KsF Brae | 
12XF6e1e//) | 

410 IF (NR-NOR)7+112,112 é ! 

4112 CONTINUE | 
END 
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APPENDIX III 


GRAPHS FOR APPROXIMATE EQUATIONS 


Graphs for the approximate equations 
used in the experimental method and 
equations section are given on the 


following pages. 
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Figure 10. 
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